Phonon Density of States of LaFeAsOi 
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We have studied the phonon density of states (PDOS) in LaFeAsOi_a:Fa; with inelastic neutron 
scattering methods. The PDOS of the parent compound (x=0) is very similar to the PDOS of 
samples optimally doped with fluorine to achieve the maximum Tc (x~ 0.1). Good agreement is 
found between the experimental PDOS and first-principle calculations with the exception of a small 
difference in Fe mode frequencies. The PDOS reported here is not consistent with conventional 
electron-phonon mediated superconductivity. 

PACS numbers: 74.70.-b, 74.72.-h, 78.70.Nx 63.20.kd 



The discovery of superconductivity at 28 K in fluo- 
rine doped LaFeAsO by Kamihara et al. and subse- 
quent reports of TcS in excess of 50 K in the broader 
RFeAsO (R=La,Ce,Pr,Nd,Sm,and Gd) family has gen- 
erated great interest in the condensed matter physics 
community [E i, i, i, % % 0, i, Although the su- 
perconductivity in the RFeAsO family is widely believed 
to be unconventional, there is no general consensus as 
to the precise mechanism behind superconductivity (e. 5. 
Ref. [lO| and references therein). Even the symmetry of 
the superconducting gap is controversial with reports of 
both one gap 11 ll , and two gaps 12, Ujl, |lj] and reports 
of line nodes [l3| as well as a fully gapped Fermi surface 

The critical temperatures of the RFeAsO family are 
exceeded only by the high-Tc cuprates, so it is nat- 
ural to ask if both materials share common physics. 
In the cuprates, while significant changes are observed 
in the phonon spectrum as a function of doping and 
temperature 16, 17, 18] the consensus is that supercon- 
ductivity does not arise solely from electron-phonon cou- 
pling. However, electron-phonon coupling docs play a 
significant role in the formation of competing charge- 
ordered ground states that appear in proximity to 
superconductivity[19[. To date there have been only lim- 
ited studies of the phonon modes in the RFeAsO family 
[20, I21I, I24I . Therefore, to provide further insight into the 
phonon behavior in RFeAsO, we have performed inelas- 
tic neutron scattering measurements on samples indepen- 
dently prepared by two groups of the parent compound 
(LaFeAsO) as well as samples optimally doped with flu- 
orine to achieve the maximum Tc (LaFeAsOo.gFo.i). We 
find little difference between the phonon frequencies of 
the parent compound and the superconducting counter- 
part. The experimental phonon density of states (PDOS) 
is in excellent agrement with the first-principle calcula- 
tions of Singh et aZ.[6t|. The main exception being the 
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FIG. 1: Neutron powder diffraction data: (a) LaFeAsO (SI) 
and (b) LaFeAsOo.sgFo.ii (S2). In both panels the first three 
lines are the sample phase, aluminum sample can, and FeAs. 
For (a) the fourth line is La203 while for (b) it is LaOF. 



observation of several Fe modes at '^10% lower energy 
in the experimental PDOS than predicted by calcula- 
tions. The experimental PDOS determined here is not 
consistent with conventional electron-phonon mediated 
superconductivity in fluorine-doped LaFeAsO. 

Members of the RFeAsO family crystallize in the Zr- 
CuSiAs structure type (space group PA/nmm). LaFeAsO 
has been found to undergo a structural phase transition 
at ^160 K to a low temperature orthorhombic state fol- 



TABLE I: Sample composition determined from Reitveld 
refinement of neutron powder diffraction data. The sam- 
ples are denoted as follows: S1(S3) LaFeAsO and S2(S4) 
LaFeAsO0.89F0.11 (LaFeAsOo.gFo.i). Phase fractions are re- 
ported as weight percent. The La203(LaOF) phases were 
found only in the parent (superconducting) samples. 



SI 



S2 



S3 



S4 



Sample 

FeAs 

La203/(LaOF) 



87.78 
6.17 
6.05 



95.17 
2.96 

(1.87) 



82.80 
8.03 
9.17 



93.10 
4.49 
(2.42) 



lowed by transition to an antiferromagnetic ordered state 
0)l23>l2J|- Doping with fluorine suppresses the structural 
distortion and transition to long range magnetic order 
and superconductivity emerges with a maximum Tc '^28 
K in LaFeAsO0.89F0.11 0, 



2l|,l24|. 



We have investigated two sets of samples which were 
synthesized independently at Oak Ridge National Labo- 
ratory and at Ames Laboratory. The Oak Ridge samples 
were synthesized by pelletizing a finely ground mixture of 
FeAs, La, La203, and LaFa (for the fluorine doped sam- 
ple) in a glove box. The resulting pellet was then sealed 
in a silica tube with 1/3 atm. argon and heated twice 
at 1200 °C for 15 h. All chemicals used were from Alfa 
Aesar and were of 99.95 % purity or better. The synthe- 
sis procedure resulted in samples of LaFeAsO (SI) and 
LaFeAsO0.89F0.11 (82) of 4.6 g each. A T^ of 27 K for 
S2 was determined by magnetic susceptibility measure- 
ments. The Ames samples were synthesized as follows: 
Polycrystalline LaFeAsO and LaFeAsOo.gFo.i samples 
were synthesized by a standard solid reaction method. 
LaFeAsO samples were prepared using LaAs, Fe2 03 and 
Fe, while the 10% F doped samples used LaAs, LaFs, 
Fe203, Fe and As. A homogeneous mixture of starting 
materials was pressed into pellets and placed in a quartz 
tube, which was then sealed under Ar partial pressure. 
The entire preparation was done in N2 filled glove boxes. 
The ampoule was kept at 600° C for 5 hr before heating to 
950°C with ramp rate 50°C/hr, held at temperature for 
15 hr, and then annealed at 1150°C for 50 hr. Finally, the 
ampoule was furnace-cooled. The synthesis process re- 
sulted in samples of LaFeAsO (S3) and LaFeAsO0.89F0.11 
(S4) of 5.0 g each. The T^ of 23 K for S4 was determined 
by magnetic susceptibility measurements. 

For single-phonon scattering in a polyatomic material, 
the intensity of inelastic scattering of neutrons from a 
polycrystalline sample has the form, 

5(g,c.) = 5]cT.^exp(-2W^0^^[«M + l], (1) 

i 

where Ui and m.i are the neutron scattering cross sec- 
tion and mass of atom i, exp(— 2Wi) is the Debye- Waller 
factor, n{ijj) -|- 1 is the Bose occupation factor. Gi{uj) is 




FIG. 2: Inelastic neutron scattering intensity (color scale) of 
LaFeAsO0.89F0.11 (S2) at 10 K as a function of Q and E. The 
scattering due to the empty sample holder has been removed. 



defined as 



G»(^) = :T^^|e.(j,q)p,5[w-w(j,q)]. 



(2) 



jq 



where w(j, q) and ei(j, q) are the frequencies and eigen- 
vectors. Equations [T] and [5] show that inelastic neutron 
scattering measures the weighted PDOS. In particular, 
the partial PDOS (Gj(w)) from each atomic species is 
weighted by ai/rrii. For the elements in LaFeAsO, the 
relative magnitudes of this quantity are 0.23, 0.76, 0.27 
and 1 for La, Fe, As, and O respectively and, hence, 
the vast majority of the scattered intensity is from Fe 
and O. This weighting makes inelastic neutron scatter- 
ing very sensitive to the presence of certain atoms and, 
consequently, understanding the levels of impurities are 
crucial. To this end, we have performed neutron diffrac- 
tion measurements at room temperature on all samples 
using the HB3 triple-axis spectrometer at the High Flux 
Isotope Reactor (HFIR) at Oak Ridge National Labora- 
tory. The results of these measurements are shown in 
Fig. 1 for SI (a) and S2 (b). Reitveld refinement of 
the diffraction data[25J indicates that in addition to the 
main phase, the samples contain small amounts of other 
phases as shown in Table lUl A La(0H)3 phase may also 
form in LaFeAsO samples, due to exposure in air. This 
phase is of concern for the PDOS measurements as the 
a /m weighting makes this technique very sensitive to hy- 
drogen in materials. We only found the La(0II)3 phase 
in S3 in very small quantities. 

The inelastic neutron scattering measurements re- 
ported here were performed on the ARCS direct geome- 
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FIG. 3: Phonon density of states of LaFeAsOi-^F^ (E,=130 
meV). (a) Experimental PDOS at 35 and 300 K for 
LaFeAsOo.gFo.i (S4). (b) Experimental PDOS for LaFeAsO 
(S3) and LaFeAsOo.gFo.i (S4) at T=300 K. (c) First- 
principles calculations of Singh et al.^. The partial density 
of states have been normalized as described in the text. 



try time-of-flight chopper spectrometer at the Spallation 
Neutron Source at Oak Ridge National Laboratory. A 
closed cycle refrigerator provided sample temperatures 
from 10 to 300 K. The data presented were measured 
with three different incident energies (EJ / chopper fre- 
quencies, 130 meV / 600 Hz, 100 meV / 600 Hz, and 
30 meV / 300 Hz. These three settings result in an en- 
ergy resolution of approximately 6.9 meV, 3.9 meV, and 
1.2 meV respectively at the elastic position. Limited in- 
elastic scans at constant Q were performed on the HB3 
triple-axis spectrometer at the HFIR (not shown) that 
provided independent confirmation of the ARCS data. 

The inelastic neutron scattering intensity of 
LaFeAsO0.89F0.il (S2) as a function of Q and E at 
10 K with an incident energy of 100 meV is shown in 
Fig. [21 The scattering from the empty holder has been 
subtracted-no other corrections have been made to the 
data. A number of excitations are clearly identifiable 
as lines of intensity at constant energy transfers with 



the strongest peaks observed near 12, 25, 30, and 40 
meV. The measured intensity becomes much stronger 
as Q increases allowing these modes to be identified as 
phonons due to the Q^ dependence of the scattering 
law (Eq. [T]). The intense line at zero energy transfer 
is the elastic line and in some cases the dispersion of 
acoustic modes is observed near Bragg points. The 
extra intensity observed at low Q is due to low angle 
instrumental background. 

Figure ^a) shows a comparison of LaFeAsOo.gFo.i 
(S4) at 35 and 300 K. The data were analyzed in the inco- 
herent scattering approximation. This analysis includes 
corrections for Bose and Debye- Waller factors, detector 
efhciency, as well as multiphonon scattering (including 
some account of multiple scattering) , as described in [20| . 
The resulting PDOS is normalized to unity in all cases. 
In Figs. [3] and |4] the dashed line indicates where the 
measured data ends and a quadratic approximation to 
the low energy PDOS begins. Fig. [Htb) shows the PDOS 
at 300 K for the parent compound (S3) together with the 
superconductor (S4) for E^ of 130 meV at 300 K. There is 
no discernable difference between the inelastic response 
of the parent and that of the superconducting sample. 
This is in stark contrast to the high Tc cuprates where 
changes in oxygen concentration in YBa2Cu307_5 leads 
to large changes in the PDOS at higher energies idllq. 

In addition to the similarity between the parent and 
superconducting samples the phonon modes in the su- 
perconducting samples show only small variation with 
temperature. The data of Fig. [Hfa) at 35 K and 300 K 
are of high statistical quality, and are similar at nearly 
all energy transfers with both data sets clearly showing 
a phonon cutoff of 68 meV. The major exception to this 
is a 1-2 meV stiffening of the 31 meV peak while cooling 
from 300 K to 35 K. This stiffening is also shown in Fig. 
|4fa) which shows higher resolution measurements for the 
superconducting sample (S2) at temperatures of 300 K, 
40 K, and 10 K. Additionally, the peak at 12 meV softens 
slightly on cooling from 300 K to 40 K. As shown in Fig. 
m no large changes are observed in cooling below Tc. 

Figurepljc) shows the results of first-principles calcu- 
lations [Jof the bare PDOS of LaFeAsO. To allow for 
more direct comparison with measurements, the partial 
calculated PDOS is weighted by the relevant a/m ratio 
and convolved with an approximation of the instrumental 
resolution for E^ of 130 meV. In general, the calculation 
agrees very well with the measured PDOS. The major 
deviation between experiment and theory is a number 
of the Fe modes which appear at systematically lower 
energies than the theoretical prediction. This is most 
evident near 30 mcV where calculations predicts an Fe 
mode near 34 meV while experiment shows a strong peak 
at 31 meV. This represents a 10% shift in the Fe mode 
energy. To further emphasize this point, the Fig. Iljb) 
shows high resolution measurements taken with E^ of 30 
meV and the results of calculations treated in the man- 
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FIG. 4: The Phonon Density of States of LaFeAsOo.ggFo.ii 
(S2). (a)Shows data with Ei=100 meV at 10, 40, and 300 K. 
(b) Shows a comparison of high resolution data collected with 
Ei=30 meV at 300 K. 



ner described above. The calculation predicts a mode at 
17.2 meV which is not observed experimentally, but as 
in the case of the Fe mode at 31 meV a 10% shift of the 
Fe mode frequencies in the calculation would explain the 
relatively larger density of states observed in the experi- 
mental PDOS near 15 meV. This result is consistent with 
an inaccurately small Fe-As interatomic distance from 
local density approximation (LDA) calculations 10|, |27| , 
resulting in the prediction of higher phonon frequencies. 
It is important to note that despite these deviations 
between experiment and theory, a 10% shift in Fe mode 
frequencies would not appreciably change the Migdal- 
Eliashberg theory calculation of Boeri et. al [28|. Al- 
ternatively, assuming that the experimental PDOS deter- 
mined with inelastic neutron scattering approximates the 
electron- pho non spectral density, a'^F(uj), (see for exam- 
ple Ref. [29|) similar values are obtained to those found 
by Boeri et al. [28| for the phonon-electron couphng 
constant, logarithmic average of phonon frequency, and 
Tc (~0.5 K). As such, these measurements do not sup- 
port conventional electron-phonon mediated supercon- 
ductivity. We note that this analysis does not take into 
account non-standard electron phonon coupling mecha- 



nisms, such as the Fermi surface nesting effects known to 
be present in these compounds [30|. 

In conclusion, we have measured the phonon density 
of states of the parent compound LaFeAsO and super- 
conducting LaFeAsOo.gFo.i. The measured phonon fre- 
quencies agrees well with first-principle calculations with 
the exception of a systematic shift of about 10% of Fe 
mode frequencies. Moreover these measurements are not 
consistent with conventional phonon-mediated supercon- 
ductivity. 
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